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Starting from the Lifshitz formula for the Casimir force between parallel plates we 
calculate the difference between the forces at two different settings, one in which 
the temperature is Ti = 350 K, the other when T2 = 300 K. As material we 
choose gold, and make use of the Drudc dispersion relation. Our results, which are 
shown graphically, should be directly comparable to experiment. As an analogous 
calculation based upon the plasma dispersion relation leads to a different theoret- 
ical force difference, an experiment of this kind would be a decisive test. We also 
present an analogous calculation for the case when the two plates are replaced with 
a sphere-plate system, still with gold as material in both bodies. The sphere is 
assumed so large that the proximity theorem holds. Discussion of the consistency 
with the third law of thermodynamics and the validity of the surface impedance 
approach is provided. 



1 Introduction 

Recently, there has been something of a heated controversy concerning the 
temperature dependence of the Casimir effect between parallel metallic plates. 
The classic result for thepressure on one ideal metal plate separated from 
another by a distance a is^ (f3 = 1/T; in this paper we set h = c = fcs = 1) 
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Recently, however, it has been suggested^ that the transverse electric zero 
mode should not contribute, which if the metal is otherwise regarded as ideal 
(reflection coefficients equal unity), would lead to a presumably observable 
linear temperature correctio n to the Casimir force, and a violation of the 
third law of thermodynamics!*^! However, real metals do not possess ideal 
reflection coefficients, and so are more complicated, and there seems to be no 



1 



contradiction with the Nernst heat theorem.^ It is the purpose of this paper 
to address the observable and thermodynamical consequences on the Casimir 
pressure of using the observed permittivity of real metals, particularly gold, 
which has been used in many recent experiments. 

2 Proof that TE Zero-Mode Does Not Contribute at T = 

The crucial observation that the transverse electric zero mode does not con- 
tribute is based on the condition \im u ^o lu 2 e(ui) = 0. This follows from the 
general dispersion relation^ 

XM = -j = A duJ — / 2 T ■ \2 - 2 ) 

Here the spectral function is positive, p(u>) > 0, and satisfies the sum rule a 

du>' p(u>') = 1. (3) 

o 

The structure of this dispersion relation allows us to make the complex fre- 
quency rotation, uj — > so 

4tt J uj' 2 + C 



x(<0 = 2/ (4) 



The proof is now immediate, because 



which uses the sum rule. Since the last integral converges to one as £ — > 0, 
the desired limit is established. Of course, this behavior is consistent with the 
Drude model, and not with the plasma model. 

3 Temperature Dependent Lifshitz Force 

We now recall the Lifshitz expression for the Casimir force between two par- 
allel nonmagnetic plates, characterized by a permittivity e(u>), and separated 
by a gap a. In the notation of Ref. ^the force per unit area can be written as 

A p~ 2y B e~ 2y 



^ t = -4tE'/ y 2d y 



(6) 



l-A m e- 2 v l-B m e- 2 v 
Here, y is a dimensionless quantity, y = qa, with 

<7=\Ai + C Cm = 27rm//3, 7=^0//?, (7) 

a The normalization of p(w) plays no role in our argument. The high-frequency behavior of 
the susceptibility given by the Drude formula is accepted, for which the spectral function 
is p(w) = (2/71)7/(^2 +7 2 )- 



and kj^ is the transverse wave vector (i.e., the component of k parallel to the 
plates). Further, we have defined the squared reflection coefficients A m and 
B m by 

A m =(^)\ B m =( S -A\ (8) 
\ep + sj \s+pj 

where the Lifshitz variables s and p are given by 



s = ^e-l+p 2 , P = q/C m . (9) 

The permittivity s(i(^ m ) is a function of the imaginary Matsubara frequency 
Cm- If the medium is non-dispersive, A rn (q) = A(p) and B m (q) — B(p). The 
prime on the sum in Eq. JHJ means that the m = term is counted with half 
weight. 

The free energy F per unit area follows from the relation T T = —dF/da: 

F = ^iti E ' / y d y H 1 - A ^- 2y ) + hi - B m e-^)] . (10) 

Note that A m refers to the TM mode, B m refers to the TE mode. 

We are looking for temperature effects in the Casimir force, and shall in 
the following focus attention on the following two temperatures, both easily 
accessible in the laboratory: 

Ti= 300K, T 2 = 350K. (11) 

Specifically, we want to calculate the difference in Casimir pressure between 
the two temperatures: 

AT = ^(350 K) - ^(300 K) = |^(300K)| - |JF(350K)| (12) 

(i.e., for convenience we let AJ- mean the difference between the magnitudes). 
This quantity depends on the values of A m and B m , which in turn depend 
on which dispersion relation is adopted for s{iC, m ). In this way we calculate a 
decisive quantity which in principle is directly comparable to experiment. 

4 Dispersion Relation 

As in our preceding paper, El we chose gold as the material of which the plates 
were composed. For this metal we have access to excellent numerical data for 
s(i() (courtesy of Astrid Lambrecht and Serge Reynaud). The data are shown 
graphically in Refs-L*- 1 ^. It turns out that for low and moderate frequencies, 
at least up to about 1.5 x 10 15 rad/s (1 eV), the data are nicely reproduced 
by the Drude dispersion relation 

e(i = l + ^ ry (13) 
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where, at room temperature, the plasma frequency lu p and the relaxation 
frequency v are equal to 

w p = 9.0eV, j/ = 35meV. (14) 

Strictly speaking, one should take into account also the temperature depen- 
dence of z/, such that v — > T): 

Here, v(T) can be calculated via use of the Bloch-Gruneisen formula, as ex- 
plained in Appendix C in Ref. 

At high frequencies, ( > 2 x 10 15 rad/s, the Drude formula gives values 
for e which are too low. 

In this context it is of interest to know: What frequency region gives the 
main contribution to the Casimir force? To analyze this point, it is convenient 
to go back to the expression @) , from which it is seen that the most important 
region is when y is of order unity, y *~ 1. Assuming that the transverse wave 
vector kj^ does not dominate in the expression Q for q, we thus get the 
condition 2irma/{3 ~ 1, or 

(16) 

If a = 1 /im, we have aT — 0.13 at room temperature, resulting in to = 1 
as the dominant frequency mode. If a = 0.5 fim, we expect a somewhat 
smeared-out distribution over the lowest integer values for m. If a = 3 /im 
or higher, the m = mode should be highly dominant: the problem becomes 
a high-temperature problem. We have done explicit numerical calculations, 
reproduced here in Table 1 for conveniencej^which confirm these expectations 
in detail. The numbers in the table are the relative percentage of each mode 
to, i.e., the quantities (-?"^/^ rT ) X 100, where denotes the mth mode 
contribution to the force: 

oo 

T T = E ( 17 ) 

m=0 

The numbers in the table are calculated from the empirical results for e(z£ m ). 

The Matsubara frequencies, in view of Eq. J7J, become 

_ J 2.47m x 10 14 rad/s for T = 300K, 
Cm _ \ 2.88m x 10 14 rad/s for T = 350 K. ^ ' 

From this it follows that from to = 1 up to about m — 6, the frequencies stay 
so low that the Drude formula can be used with confidence. And from Table 1 
we see that in particular for large distances, from a = 1 /im and upwards, these 
frequencies encompass the large majority of the mode contributions to the 
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Table 1. Contribution from the various Matsubara frequencies for gold. The relative con- 
tribution of T T is given, in percent, for each mode in the interval m g [0,7]. Room 
temperature is assumed. Empirical permittivities inserted, for all frequencies. 
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force. For small gap widths, a = 0.5 /xm and lower, the situation may be more 
questionable as indicated by the first line in the table, but it seems that even 
in these cases we can use the Drude formula with sufficient accuracy to make 
a meaningful comparison with experiments, given the present experimental 
accuracy. In the following we will use the Drude formula throughout. This 
simplifies the calculation significantly. 

5 Calculated Results 

Fig-Hshows how the Casimir force itself varies with a, for parallel gold plates. 
Since the force according to Eq. © is negative for attraction, we choose to 
present the magnitude \T T \ in the figure, which for convenience is presented 
in semilog form. The curve is calculated for T = 300 K, but a similar curve 
calculated for T = 350 K turns out to be visually indistinguishable from the 
curve shown. 

Fig. El shows the force difference A J 7 , calculated according to Eq. Ijl2|l ■ 
Typical orders of magnitude for AT are in the millipascal range, when a is less 
than about 0.5 fim. It is notable that it is the low-temperature term, T = 300 
K, which yields the strongest force at low and intermediate distances. (When 
a > 2.8 /j,m, AT changes sign.) We note that the influence from temperature 
in the formula © is rather complex since T occurs at three different places: 
(i) in the prefactor; (ii) in the lower limit of the integral, and (iii) in the 
dependence of A m and B m on T via the permittivity e = e(i2irmT). (We 
will discuss the temperature dependence of the relaxation frequence v below.) 
Actually we showed the essentials of the temperature dependence of the force 
already in Fig. 5 in Ref. ^1 Assuming the Drude formula, and ignoring the 
temperature dependence of the relaxation frequency, we found the force \T T \ 
to diminish with increasing T up to aT ~ 0.35. For higher values of aT, 
the force was found to increase again. In view of this, indeed we expect 
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Figure 1. Magnitude of Casimir pressure between parallel plates, versus gap width a, when 
T = 300 K. The curve for T = 350 K overlaps the one shown. 



that the lower-temperature term dominates in Fig. [21 for low and intermediate 
distances. (Note also that the change in sign for a — 2.8 /im, mentioned 
above, corresponds nicely to the expected transitional region since in this 
case aT ~ 0.37.) 

So far, we have assumed a constant relaxation frequency, v — 35 meV. 
Will our results be changed significantly if we take into account the tempera- 
ture dependence of v\ cf. Eq. (|15fl ? The answer turns out to be no. We have 
made an explicit calculation of this, based upon the Bloch-Griineisen formula, 
yielding v = 35.6 meV for T = 300 K, and v = 41.8 meV for T = 350 K. 
The results were visually indistinguishable from those given in Fig. [2 so that 
our conclusion is that the assumed constancy of v is justified for practical 
purposes. 

Fig. [3] shows how the corresponding difference in free energy, 

AF = F(350K) - F(300K) = |F(300K)| - |F(35QK)|, (19) 

varies versus a for parallel plates. As expected, it is again the case T = 300 K 
which is the dominant one for small and moderate gap widths. 

Finally, we consider the case which is probably the one of principal exper- 
imental interest, namely a sphere of radius R situated above a plane surface. 
Let the minimum distance between the spherical surface and the plane be a. 
Assuming the sphere large enough for the proximity theorem to be validP we 
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Figure 2. Force difference between parallel plates, Eq. 1121 . versus gap width. 
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Figure 3. Free energy difference for parallel plates, Eq. 1191 . versus gap width. 



have for the force J^ s on the sphere: 



= 2nRF(a), 



(20) 
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Figure 4. Force difference between sphere and plate, Eq. 1211 . versus gap width a. 



where F(a) is the free energy for parallel plates as given by Eq. (|10|l . Fig. 21 

shows how the difference 
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D A^ ps = - [^ pfl (350K) - ^ ps (300K)] 
H K 



(21) 



varies with a. Again, it is the lower-temperature term that is the dominant 
one. 

A dedicated experiment to look for the temperature dependence we have 
proposed is probably essential to settle this issue. The recent experiment by 
Decca et al.^Ms claimed to be in disagreement with our prediction. However, 
that comparison is in fact not based on our detailed calculations, and the 
experiment is subject to large, uncontrolled errors!^ 



6 Behavior of the Free Energy at Low Temperature 

The low temperature correction is dominated by low frequencies, 6 where the 
Drudc formula is extremely accurate. Using this fact, we have performed 
analytic and numerical calculations which show the free energy has a quadratic 
low-temperature dependence, independent of the plate separation: 

LU 2 

F(T) = F +T 2 — 2-(21n2- 1) = F + T 2 (19 eV), (22) 

b This statement is in the context of using of the Euler-Maclaurin summation formula to 
evaluate Eq. J§J, for example. 
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Figure 5. The behavior of the free energy for low frequencies, in the Drude model, 
with parameters suitable for gold, and a plate separation of a = 1 fim. Here _F TE = 



putting in the numbers for gold, rather than the naive extrapolation 



F = F a + T- 



C(3) 



T 



r0.30 



(23) 



167ra 2 Ana 2 

We see from Fig. [S] that this value indeed results if one extrapolates the 
approximately linear curve there for £et > 0.25 to zero, following the argument 
given in Eq. (2.8) of Ref. However, we see that the free energy smoothly 
changes to the quadratic behavior exhibited in Eq. (|22|) . 

Results consistent with these have been reported by Sernelius and 
Bostrom.EH 



7 Surface impedance form of reflection coefficient 

It has been proposed that the resolution to the temperature problem for the 
Casimir effect is that the surface impedance form of the reflection coefficients 
should be used in the Lifshitz formulapl rather than that based on the bulk 
permittivity. Here we show that the two approaches are in fact equivalent, 
and that the former must include transverse momentum dependence. 
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For the TE modes, the reflection coefficient is given b}^ 



r TE = _ h z fc 2 ^ (24) 



where 



k\ - zVC 2 [e(zC)-l] + 9 2 , (25) 



with q 2 = k\ + ( 2 , and the subscripts 1 and 2 refer to the metal and the 
vacuum regions, respectively. Now from Maxwell's equations outside sources 
we easily derive just inside the metal (the tangential components of E and B 
are continuous across the interface) 

- ik u k± • B ± - iue (l - ^] k r (nx Ej_) = 0, (26) 



-ik lz k± ■ (n x E_l) - iu;k± - B± = 0. (27) 

Here n is the normal to the interface. Now the surface impedance is defined 
by 

Ej_ = Z(oj, k x )B ± x n, n x E ± = Z(uj, k ± )B ± . (28) 
So eliminating B x using this definition we find two equations: 

*i* = -|, ( 29 ) 
k 2 lz =uj 2 e-ki, (30) 

the latter being the expected dispersion relation l|25|l . Substituting this into 
the expression for the reflection coefficient we find 



r 



TE 



C + Zq 1 + Zp q 



(31) 



C-Zq 1-Zp' C 

which apart from (relative) signs (presumablyjust a different con vention 
choice) coincides with that given in Geyer et al.El or Bezerra et alPl 



7.1 Dependence on Transverse Momentum 

However, it is crucial to note that the "surface impedance" so defined depends 
on the transverse momentum, 



VC 2 [e«)-l]+9 2 ' 

and so r TE — > as £ — * just as in the dielectric constant formulation. Of 
course, we have exactly the same result for the energy as before, since this is 
nothing but a slight change of notation. 
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It is therefore incorrect to assume that Z is only a function of frequency, 
not of transverse momentum, and to use the normal and anomalous skin effect 
formulas derived for real waves impinging on imperfect conductors. c 

How does the usual argument go? The normal component of the wavevec- 
tor in a conductor is given by 

1/2 

— > \J i4iruja, (33) 

from which the usual normal skin effect formula follows immediately, 

Z(.) = -(l-i) ] [£. (34) 

However, the last step here consists in omitting two "small" terms: e (okay) 
and k\ < ui 2 . Here this last is not valid because in going to finite temper- 
ature we have severed the connection between uj — > iC, and the latter 
is in no sense ignorable as we take £ — > to determine the low temperature 
dependence. This is the same error to which we refer in our publishe d pa per.^ 
These considerations are consistent with those of Esquivel et alM^ 



e + 1- 



.47TC7 



8 Conclusions 

Our results of main interest are probably those shown in Figs. [21 and These 
force curves show a dependence upon a that reflect our underlying choice of 
the Drude dispersion rel atio n. We can compare our results with those recently 
obtained by Chen et alPD They make use of the plasma dispersion relation 
instead of the Drude relation, and obtain results for the Casimir forces that 
differ from ours even in sign. The force curves thus give rise to a very useful 
critical test, in principle. It would be quite interesting if the experimentalists 
could measure these force curves directly. We have also commented on the 
purported violation of basic principles of thermodynamics and on the claimed 
necessity to use surface impedances at ~ 1 /im plate separations and find no 
merit in these objections. 
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c Of course, in general, the permittivity will be a function both of the frequency and the 
transverse momentum, e(tj,kx), but we believe the latter dependence is not significant for 
separations larger that hc/u) v = 0.02 fim. 
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